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Abstract 



The scalar and tensor polarizabilities of the deuteron are calculated using 
the recently developed effective field theory that describes nucleon-nucleon 
interactions. Leading and next-to-leading order contributions in the pertur- 
bative expansion predict a scalar electric polarizability of (%eo = 0.595 fm 3 . 
The tensor electric polarizability receives contributions starting at next-to- 
leading order from the exchange of a single potential pion and is found 
to be a£2 = —0.062 fm 3 . We compute the leading contributions to the 
scalar and tensor magnetic polarizabilities, finding (3mo = 0.067 fm 3 and 
(3 M 2 = 0.195 fm 3 . 
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I. INTRODUCTION 



Efforts to develop a systematic treatment of nucleon-nucleon interactions PHH1 have 



culminated in an effective field theory with consistent power counting ||19|| . The leading 
contribution to two nucleons scattering in an S-wave comes from local four-nucleon operators. 
Contributions from pion exchanges and from higher derivative operators are suppressed by 
additional powers of the external nucleon momentum and by powers of the light quark 
masses. The technique successfully describes the NN scattering phase shifts up to center- 
of-mass momenta of p ~ 300 MeV per nucleon (T^| in all partial waves. 

To accommodate the unnaturally large scattering lengths in S-wave nucleon-nucleon 
scattering, fine-tuning is required, which in turn can complicate power counting in the 
effective field theory. Dimensional regularization with power divergence subtraction (PDS), 
described in |19| , provides a consistent power counting scheme. Since the deuteron is the 
lightest nucleus and does not have irreducible forces between three or more nucleons, it 
provides a unique laboratory for studying the strong interactions. Being bound by only 
2.2 MeV, the characteristic momentum of the nucleons in the deuteron is ~ 40 MeV and 
should be well described by the effective field theory which is valid below the scale Ann ~ 
300 MeV |19| . The electromagnetic moments and form factors of the deuteron have been 



explored with this new effective field theory ^{J . The charge radius and form factor receive 
contributions from leading and next-to-leading (NLO) orders in the expansion with the 
theoretical result reproducing the measured values within the uncertainty coming from the 
omission of higher order terms. The magnetic moment and form factor receive contributions 
at leading order from only the nucleon magnetic moments. At next-to-leading order the 
deuteron magnetic moment determines a combination of counterterms that appear in the 
Lagrange density at this order. The quadrupole moment first appears from the exchange 
of a single potential pion. Pion-pole contributions, multiple potential pion exchange, and 
higher dimension operators contribute only at higher orders in the expansion. 

Unlike the electromagnetic form factors, the electric and magnetic polarizabilities of an 
object are a direct measure of its "deformation" due to the presence of external electric and 
magnetic fields. Extensive experimental and theoretical progress has been made in under- 
standing the electric and magnetic polarizabilities of the nucleon (for an overview see 
Chiral perturbation theory provides a systematic theoretical analysis in the single nucleon 



sector, with one- loop pion graphs dominating the electric polarizability, e.g. p2|-[26|1. The 
magnetic susceptibility, on the other hand, is dominated by the A-pole and has a signifi- 
cant uncertainty associated with it, e.g. [^7J]. Recently, the discussion has been extended 
to include "generalized polarizabilities," the amplitudes appropriate for interactions with 
electrons 



Theoretical understanding of the polarizability of the deuteron has been expressed in 
terms of meson exchange potential models |3~T^f47}| (for an excellent discussion see JJ7|). In 
this work we present a model independent, analytic computation of the electric polarizabili- 
ties of the deuteron to NLO and the magnetic polarizabilities to leading order in the effective 
field theory describing nucleon-nucleon interactions. 
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II. EFFECTIVE FIELD THEORY FOR NUCLEON-NUCLEON INTERACTIONS 



The terms in the effective Lagrange density describing the interactions between nucleons, 
pions, and photons can be classified by the number of nucleon fields that appear. It is 
convenient to write 

£ = jC + jC 1 + jC 2 + (2.1) 

where C n contains n-body nucleon operators. 

Cq is constructed from the photon field A^ 1 = (A , A) and the pion fields which are 
incorporated into an £77(2) matrix, 

where / = 132 MeV is the pion decay constant. S transforms under the global SU(2) L x 
SU(2) R chiral and U(l) em gauge symmetries as 

S -> LEi^, S -> e ia0 -"Ee _ * QQem , (2.3) 

where L G SU(2)l, R G SU(2)r and Q errt is the charge matrix, 

g — = ( ) ' (2 - 4) 

The part of the Lagrange density without nucleon fields is 

1 f 2 f 2 

C = -(E 2 -B 2 ) + J —Tr D^ZD^ + ±-ATr m g (E + E f ) + ... . (2.5) 

2 o 4 

The ellipsis denote operators with more covariant derivatives D^, insertions of the quark mass 
matrix m q = diag(m u , m^), or factors of the electric and magnetic fields. The parameter A 
has dimensions of mass and m 2 = \{m u + m d ). Acting on S, the covariant derivative is 

DfjT, — d^Y, + ie[Q em , . (2.6) 

When describing pion-nucleon interactions, it is convenient to introduce the field £ = 
exp (ill//) = \/E. Under SU(2) L x SU{2) R this transformations as 

£ L£[/t = ^ i2 t ) (2.7) 

where U is a complicated nonlinear function of L, R, and the pion fields. Since U depends on 
the pion fields it has spacetime dependence. The nucleon fields are introduced in a doublet 
of spin 1/2 fields 

N = (2 ' 8) 

that transforms under the chiral SU(2) L x SU(2) R symmetry as — > UN and under the 
U(l) em gauge transformation as — * e iaQem N. Acting on nucleon fields, the covariant 
derivative is 
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D^N = (<9 M + V^ + ieQ em A„)N 



(2.9) 



where 



(2.10) 



The covariant derivative of N transforms in the same way as N under SU(2)i x SU(2)r 
transformations (i.e. D^N — > UD^N) and under U(l) gauge transformations (i.e. D^N — > 



e iaQem D^N). 



The one-body terms in the Lagrange density are 



& = NU iD + 



D 2 



2M 



2M 



N , 
«1 



iV - y^vV ■ (£D£ f - £ f D£)iV 



+ 2 7 r4 Aro) iVtiVE 2 + 2vr4 Ari) iVtr 3 iVE 2 + 2tt (3^ 0) NB 2 + 2-k^ 1 ^WB 2 + . . . , (2.11) 

where n = ^(k p + K n ) and isoscalar and isovector nucleon magnetic 

moments in nuclear magnetons, with 



k p = 2.79285 



-1.91304 



(2.12) 



The isoscalar and isovector electric polarizabilities of the nucleon are a^°^ and a^ v±; while 
the corresponding magnetic quantities are (3^°^ and Pm • Experimentally, it is found that 



a p = (12.1 ± 0.8 ± 0.5) x 10~ 4 fm 3 



ftv 



(2.1 =)= 0.8 =F 0.5) x 10~ 4 fm 3 



(2.13) 



for the proton |}41] , [42| while the two measurements of the neutron electric polarizability 
a n = 12.0 ± 1.5 ± 2.0 [^3] and a n = ± 5 [H indicate a sizable uncertainty. 
The two-body Lagrange density is 



C { o Sl) + D^°A'L:mg))(JV T P i iV)t(iV T P i JV) 



+ 



N + h.c. 



(N T PiNy \^N T PiT) +5 Pi- 2t5PiI) 

+eL 1 (A^ t (j • BN){N^N) + eL 2 {N^a ■ Br a N){N^T a N) 
+ 2na 4 (N T PiN)^ (N T PiN)E 2 + 2vr/3 4 (iV T P i iV) t (iV T P i iV)B 2 + 



(2.14) 



where Pi is the spin-isospin projector for the spin-triplet channel appropriate for the deuteron 



P = 



1 



(2.15) 



The a matrices act on the nucleon spin indices, while the r matrices act on isospin indices. 
The local operators responsible for S — D mixing do not contribute at either leading or NLO. 
The Cf 5l) , Cf 5l) and L>f Sl) coefficients are determined from NN scattering to be 
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cf Sl) (m^) = -5.51 fm 2 , D { * Sx \m v ) = 1.32 fm 4 , c£ Sl \m w ) = 9.91 fm 4 , (2.16) 



where we have chosen to renormalize the theory at a scale fi = in the PDS scheme |L9|]. A 
linear combination of the coefficients L\p contribute to the magnetic moment of the deuteron, 
but neither they nor the coefficients «4 and (3^ contribute to the deuteron polarizability at 
the order to which we are working. The four individual nucleon polarizability counterterms 
terms also do not contribute at the order to which we are working. In eq. ( p. 14] ) we have 



only shown the leading terms of the expansion in meson fields, namely the terms we need 
for our leading plus NLO calculation. 

Since we are working with a field theory, no ambiguities arise from how we choose to 
define the interpolating fields. We can consistently neglect all operators that vanish by the 
equations of motion |45|J46]| . An operator that vanishes by the equations of motion makes a 
contribution to an observable that has exactly the same form as higher dimension operators 
that are present in the theory. Further, such operators can be removed by field redefinitions, 
and therefore it is consistent to work with a Lagrange density that does not contain operators 
that vanish by the equations of motion. 

Another point of interest is that we do not need to include the A as an explicit degree of 
freedom in the theory. In order to have a theory that is well defined for processes involving 
momenta up to ~ 1 GeV, the A must be included as a dynamical object [ffTfl . However, 
the power counting for the effective field theory describing the nucleon-nucleon interaction 
outlined in this section is limited to momenta less than p ~ A^v ~ 300 MeV. The momen- 
tum scale making the dominant contribution to graphs involving the A is approximately 



p ~ y Mn(Ma — Mn) ~ 500 MeV, higher than K^n- Therefore, the A is not included as 
a dynamical object, but its effects are included in the coefficients of the local operators (as 
are the effects of all particles not included as dynamical degrees of freedom). 



III. COMPUTING THE POLARIZABILITIES OF THE DEUTERON 



The Lagrange density described in the previous section in terms of nucleon field operators 
matches onto an effective theory describing the dynamics of the deuteron field 



C D =V\(iD Q + 



D 5 



V a - ifi D e abc VlV b B 



2M DJ 

+ 2na E0 V\V a E 2 + 2tt/3 M o V\V a B 2 



+ 2ixa E 2 
+ 2vr/3 M2 



V\V b + V\V a - -5 ab V\W 



V\V b + V\V a - ~5 ab V]p< 



E a E 6 



B a B b 



(3.1) 



T> a is an operator that annihilates a deuteron, and its spin index takes values a = 1, 2, 3. The 
covariant derivative acting on the deuteron field is = + ieQA^. The coefficient /i_d is 
the deuteron magnetic moment and has been determined from eqs.( |2.1l"D and ( |2.14|) to NLO 
in [^0J. The scalar electric polarizability of the deuteron is o:eo, and the scalar magnetic 
polarizability is /3mo- These operators give rise to interactions that are independent of the 
alignment of the deuteron with respect to an applied electromagnetic field. The tensor 
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electric polarizability of the deuteron is aE2 } and the tensor magnetic polarizability is (3m2- 
These operators give rise to interactions that depend upon the alignment of the deuteron 
with respect to an applied electromagnetic field. In order to compute q:eo,E2 and /3mo,M2 we 
will use the formalism developed in [p0[| . The polarizabilities a and f3 each have perturbative 
expansions in powers of Q = tji^/Ann ~ VMnB / Ann, where B is the deuteron binding 
energy. We denote the contribution from a given order by a superscript, e.g. the scalar 
electric polarizability 

meo = a { E~o ] + ^eo + •• > (3- 2 ) 

and similarly for the other polarizabilities. 

In matching onto the deuteron effective Lagrange density in eq. ( |3.1| ) from the Lagrange 
density describing nucleon dynamics, eqs. ( p.llj ) and (|2.14|) , we will recover the coefficient 
of each operator order by order in the Q expansion, including the coefficient of the operator 
T>] ] G 2 T> a . At leading order in the loop expansion, we find this operator to have a coefficient 
4]^, while at next-to-next-to-leading order (NNLO) it will have a coefficient + g^r- 

When determined to all orders it will sum to 2 {2M n -b) = 2Tfc' wn * c h we have written in 



eq. ( |3.1| ). However, it is important to realize that we will only recover interactions coming 
from this operator order by order in the expansion. 

It is instructive to begin by power counting the contributions from operators that appear 
in the Lagrange density eqs. ( |2.11j) and Q2.14Q . Consider the single loop graph of Fig. [I] 
between a source that creates a spin triplet NN pair and one that annihilates it, minimally 
coupling to two photons. In terms of the expansion parameter Q, a non-relativistic loop 
integral scales as Q 5 , a non-relativistic nucleon propagator as Q~ 2 , and a gradient operator 
as Q l . In order to match onto the E 2 or B 2 operator in the deuteron Lagrange density 
( |3.1|) , we must find the coefficient of uj 2 or k 2 by expanding the graph in powers of uj and 
k, respectively, where (co>,k) is the photon four-momentum. Therefore, this diagram will 
scale like Q~ l J2n { Mn q2^ ) • Wave function renormalization introduces a factor of Q, 
and therefore we find that the leading contribution to the electric polarizabilities arising 
from this graph is of order Q~ 4 . This power counting also shows that the C^/i) operator 
and the exchange of a single potential pion contribute at order Q~ 3 after wave function 
renormalization. The magnetic polarizabilities receive contributions from all the graphs 
shown in Fig. [l] at order Q~ 2 . The four-nucleon polarizability counterterms, and fa, 
appearing in (|2.14| ) contribute at order Q 1 , and can be safely neglected. Counterterms for 
the individual nucleon polarizabilities appearing in (|2.11|) contribute at order Q°. Meson 
loop corrections that are the dominant contribution to the electric polarizability of the 
nucleon, for instance Fig. |2|, appear at order three orders higher in the expansion than 
the leading order contributions. Therefore, it is probable that the polarizabilities of the 
individual nucleons will not be extracted from the polarizabilities of the deuteron. 

We will compute the leading and NLO contributions to the scalar and tensor electric 
polarizability of the deuteron, and the leading contributions to the magnetic scalar and 
tensor polarizabilities. The loop graphs in Fig. [l] give the leading contribution to the deuteron 
polarizabilities, which after wave function renormalization are 

(_ 4) _ aM N 

a E0 ~ 3^4 
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FIG. 1. Leading order contributions to the deuteron polarizabilities. The crossed circles denote 
operators that create or annihilate two nucleons with the quantum numbers of the deuteron. The 
dark solid circles correspond to the photon coupling via the nucleon kinetic energy operator ( minimal 
coupling), while the light solid circles denote the nucleon magnetic moment operator. The solid 
lines are nucleons. Only the graph with the nucleons minimally coupled to the electromagnetic field 
contributes to the electric polaribility. The photon crossed graphs are not shown. The bubble chain 
arises from insertions of the four nucleon operator with coefficient Cq . 
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+ 
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6n 



> (1) ) 2 M JV ^ 0)( 



27T7 



7 



-B) 



(3.3) 



where 7 = \/MnB and a = 1/137 the fine structure constant. The scattering amplitude in 
the ^Sq channel for a centre of mass energy E is 



(E) 



-a 



^-M N E - ie) 



(3.4) 



which is of order Q . The coefficient Cq has been determined from nucleon- nucleon 
scattering in the ^0 channel to be ]19[ Cq So ^ = —3.34 fm 2 . 

At NLO there are contributions from the exchange of a single potential pion, Fig. and 
from the operator with coefficient C*2(/i), Fig. |j. The operator with coefficient D2{[i) does 
not contribute to the polarizability of the deuteron. We find that at order Q~ 3 
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FIG. 2. A contribution to the deuteron electric polarizability from a graph that also contributes 
to the polarizability of the nucleon. The crossed circles denote operators that create or annihilate 
two nucleons with the quantum numbers of the deuteron. The solid circles correspond to the photon 
coupling via the nucleon kinetic energy operator (minimal coupling). Dashed lines are pions and 
solid lines are nucleons. 



( _ 3 ) aM%„,^,„ _ A2 ag\M 2 N m 2 (3m 2 + i6m, 7 + 24 7 2 



64tt7 3 u 3847r/ 2 7 3 (m 7r + 2 7 ) 4 

( _ 3 ) ag 2 A M 2 N 2ml + H^7 + 16m^ 7 2 + 8 7 3 
am ~ 80vr/ 2 7 2 (m. + 2 7 )4 • 

The leading contribution to the tensor polarizability of the deuteron aE2 comes from the 
exchange of a potential pion. The explicit renormalization scale dependence that appears in 
eq. (|3.5| ) is compensated by the \i dependence of the coefficient C^/f), which scales ~ l//x 2 
0- 

Numerically, we find that the electric polarizabilities are 

n - J' 4 ) , ^,(-3) , 
®E0 — &E0 + a E0 + ••• 

= 0.386 fm 3 + (0.153 + 0.057) fm 3 + ... 

= 0.595 fm 3 + ... , (3.6) 

where the dots denote contributions higher order in the expansion. The first term in paren- 
thesis on the second line of eq. Q3.6Q comes from the C^/i) operator while the second term 
arises from the exchange of a single potential pion. The convergence of the expansion for a^o 



appears to be approximately the same as it is for the static electromagnetic moments |20 
with each order being suppressed by between 1/3 and 1/2, consistent with the A^tv expan- 
i9|. We expect that the uncertainty in this numerical value is roughly Aago ~ 0.1 fm 3 



sion 

from an estimate of the next order contribution. For the tensor polarizability we find 

a E 2 — a E2 + a E2 + ... 

= -0.062 fm 3 + ... , (3.7) 

where we recall a^ 4 ^ = 0. The fractional uncertainty in otE2 is much greater than that for 
oeo as it has a vanishing leading order contribution, and we naively estimate an uncertainty 
of ActE2 ~ 0.03 fm 3 . The leading contribution to the magnetic polarizabilities are 
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FIG. 3. Graphs from potential pion exchange that contribute to the deuteron polarizabilities at 
NLO. The crossed circles denote operators that create or annihilate two nucleons with the quantum 
numbers of the deuteron. The solid circles correspond to the photon coupling via the nucleon or 
meson kinetic energy operator (minimal coupling) or from the gauged axial coupling to the meson 
field. Dashed lines are mesons and solid lines are nucleons. Photon crossed graphs are not shown. 



(5 { M l ] = 0.067 fm 3 , 

d^f = 0.195 fm 3 . (3.8) 

The large values for the magnetic polarizabilities come from the isovector magnetic moment 
of the nucleon, K\. 

It is informative to decompose these interactions into the polarizabilities of the individual 
magnetic substates of the deuteron. The electric polarizability of the m — ±1 states of the 
deuteron is a_eo — §«£2 while the polarizability of the m = state is a^o + |«_B2, and 
similarly for the magnetic polarizabilities, we find 

a l ™ l=1 = 0.637 fm 3 + ... , /?j^ l=1 = -0.063 fm 3 + ... 
a™=° = 0.511 fm 3 + ... , f3^=° = 0.327 fm 3 + ... , (3.9) 

where the ellipsis denotes higher order contributions. Numerical evaluation of our analytic 
results agree within uncertainties with values for the electric polarizabilities obtained from 
potential models. To a very high degree of precision potential models predict a scalar electric 
polarizability of (Xeq — 0.6328 ± 0.0017 fm 3 [37[], which is remarkably consistent with the 
"zero-range" limiting value of a#o = 0.632 ±0.003 fm 3 ||31|| . Further, the calculations of EC 



find a ™ 1 = 0.669 fm 3 and ct™ =0 = 0.555 fm 3 . We expect to deviate from these values 
at higher orders in the effective field theory expansion as generally potential models do not 
properly describe chiral dynamics, pion propagation, and relativistic effects. The power of 
the effective field theory formalism is that such effects can be included systematically and 
the natural size of higher order terms is known. Effective field theories will never match the 
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oo 




<dboo 

oc^> 

FIG. 4. Graphs from insertions of the operator with coefficient C2{n) that contribute to the 
deuteron polarizabilities at NLO. The crossed circles denote operators that create or annihilate two 
nucleons with the quantum numbers of the deuteron. The solid circles correspond to the photon 
coupling via the nucleon kinetic energy operator (minimal coupling) while the solid square denotes 
the C^Qu) operator. The solid lines are nucleons. Photon cross graphs are not shown. 




precision of potential or other models because their reliance on only the symmetries of the 
underlying theory and experimental results necessarily makes them less restrictive. However, 
the fact that effective field theory descriptions provide a model independent treatment makes 
them an important tool in the study of low energy processes. 



IV. CONCLUSIONS 

We have performed a model independent computation of the scalar and tensor polariz- 
abilities of the deuteron in the effective theory describing nucleon-nucleon interactions. The 
scalar electric polarizability receives leading order contributions while the tensor electric 
polarizability begins at NLO in the power counting, and so is sizably smaller than the scalar 
polarizability. The tensor electric polarizability arises only from the exchange of potential 
pions at NLO. The leading contribution to the magnetic polarizability of the deuteron comes 
from the magnetic moments of the individual nucleons in addition to their minimal coupling 
to the electromagnetic field. From this analysis we conclude that it is unlikely that the 
neutron polarizabilities will be extracted from the deuteron polarizabilities. 

We are encouraged that this effective field theory calculation reproduces the results of 
potential model calculations. The benefits of using an effective theory are that closed form 
expressions have been obtained for the polarizabilities, and through the universal couplings 
in the Lagrange density we can see how QCD relates a variety of two nucleon processes in a 
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systematic way. Further, higher order terms not considered in this work are parametrically 
smaller that those we have presented. In subsequent work we will present the photon- 
deuteron scattering cross section over a range of photon energies. 



This work is supported in part by the U.S. Dept. of Energy under Grants No. DE-FG03- 
97ER4014 and DE-FG02-96ER40945. 
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